Finite-Size Bosonization of 2-Channel Kondo Model: a Bridge between Numerical 
Renormalization Group and Conformal Field Theory 
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We generalize Emery and Kivelson's (EK) bosonization-refermionization treatment of the 2-channel 
Kondo model to finite system size and on the EK-line analytically construct its exact eigenstates and 
finite-size spectrum. The latter crosses over to conformal field theory's (CFT) universal non-Fermi- 
liquid spectrum (and yields the most- relevant operators' dimensions), and further to a Fermi- liquid 
spectrum in a finite magnetic field. Our approach elucidates the relation between bosonization, scal- 
ing techniques, the numerical renormalization group (NRG) and CFT. All CFT's Green's functions 
are recovered with remarkable ease from the model's scattering states. 

PACS numbers: 72.15.Qm, 75.30. Hx, 71.10.Hf, 11.25.Hf 

defined by H = H + H z + Hj_ + H h (ft = v F 



A dynamical quantum impurity interacting with me- 
tallic electrons can cause strong correlations and some- 
times lead to non- Fermi-liquid (NFL) physics. A proto- 
typical example is the 2-channel Kondo (2CK) model, in 
which a spin- 1/2 impurity is "overscreened" by conduc- 
tion electrons, leaving a non-trivial residual spin object 
even in the strong-coupling limit. Many theoretical treat- 
ments of this model have been developed ffl, including 
Wilson's numerical renormalization group (NRG) [|]]|] 
for the crossover from the free to the NFL regime, Af- 
fleck and Ludwig's (AL) conformal field theory (CFT) 
HQ for exact thermodynamic and transport quanti- 
ties, valid only near the NFL fixed point, and Emery 
and Kivelson's (EK) bosonization-refermionization map- 
ping onto a resonant-level model [ 5| , valid on a line 
in parameter space that connects fj the free and NFL 
fixed points. In this Letter we elucidate the well-known 
yet remarkable fact that these three approaches, despite 
tremendous differences in style and technical detail, yield 
mutually consistent results: we show that EK bosoniza- 
tion in a system of finite size L yields NRG- like finite-size 
spectra, and reproduces all known CFT results. 

Our method requires no knowledge of CFT, only that 
we bosonize and refermionize with care: Firstly, we con- 
struct the boson fields <fi and Klein factors F in the 
bosonization relation ip ~ Fe~ 1 ^ explicitly in terms of 
the model's original fermion operators {ckaj}- Secondly, 
we clarify how the Klein factors for EK's refermionized 
operators act on the original Fock space. Thirdly, we 
keep track of the gluing conditions on all allowed states. 
This enables us (i) to explicitly construct the model's 
finite-size eigenstates; (ii) to analytically obtain NRG- 
like finite-size spectra that cross over from free to CFT's 
universal NFL spectra; (iii) to describe magnetic-field- 
induced crossovers exactly; and (iv) for L — > oo to easily 
recover all AL's CFT results f§. 

The model. — We consider the standard anisotropic 
2CK model with a linearized energy spectrum 
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Here c£ ■ creates a free-electron state \kaj) with spin 
a = (|, I), flavor j = (1,2) = (+,—), radial momen- 
tum k = \p\ — pp, and normalization {ck a j, c k' a 'j'} = 
Skk'S aa 'Sjj' . We let the large- 1 k \ cut-off go to infinity 
and quantize k by defining 1-D fields with, for simplicity, 
antiperiodic boundary conditions at x = ±L/2 B, 
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where k = Al (nk — 1/2) and = 27r/L is the mean 
level spacing. By : : we denote normal ordering relative 
to the Fermi ground state |0)o- H z + H± is the Kondo 
coupling (with dimcnsionlcss A 2 =/= Xj_ = X x = X y ) to 
a local spin-1/2 impurity S a (with ^-eigenstates | 
l-D-)), and Hh describes magnetic fields hi and h e coupled 
to the impurity spin and the total electron spin Af s . 

Conserved quantum numbers. — Diagonalizing H re- 
quires choosing a suitable basis. Let any (nonunique) si- 
multaneous eigenstate of N a j = ^2k :c laj Cka 3 •> counting 
the number of (aj) electrons relative to |0)o, be denoted 
by \N) = |iV n )®|A^i)®|iV T 2)<g>|A^2}, withA^eZ 4 . Since 
H conserves charge, flavor and total spin, it is natural to 
define new counting operators, Af y (y — c, s, /, x), 



/Af c \ 

Af f 
\AT X J 



f l 1 

1 -1 

1 



y 

1 1 

1 -1 
1 -1 -1 





(N n \ 




Nix 




N n 


) 


\N l2 J 



(2) 



which give the total electron number, the electron spin, 
flavor, and spin difference between channels, respectively. 
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Equation (|J) implies that the eigenvalues N are either 
all integers or all half-integers (i.e. JV € (Z + P/2) 4 , with 
P = (0, 1) for even/odd total electron number), and that 
they obey the free gluing condition 

N c ±N f = {N s ± AQ mod 2 . (3) 

All non-zero matrix elements of H± have the form 
(A/" C) Sr-f,A/},A4;ft \H±lAr c ,Sr + l,Mf,K ± l-A), 
and since the total spin St = Ms + S z is conserved, 
the A» s -eigenvalue flips only between StT\, i- e - it fluc- 
tuates only "mildly". In contrast, the M x -eigenvalue 
fluctuates "wildly" [an appropriate succession of spin 
flips can produce any Af x satisfying (fj)]; this will 
be seen below to be at the root of the 2CK model's 
NFL behavior (in revealing contrast to the 1CK model, 
which has no wildly fluctuating quantum number, 
and lacks NFL behavior). For given (Af c , St, A/"/) 
it thus suffices to solve the problem in the correspond- 
ing invariant subspace J^e'Af |A/" C , St — \ , A// , Af x ; ft) 
|A4, St + 5, A//, A4+ 1; JJ.) , to be denoted by S, where the 
prime on the sum indicates its restriction to A/^-values 
respecting (||). 

Bosonization. — To bosonize Q the model in terms of 
the original Ckaj 's |8|,[§ , we define bosonic fields through 
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which account for particle-hole excitations (the 6's by 
construction satisfy [b qa j, b q / a ,j,] = S qq )5 aa >o~jj> and 
[b qa j,N a iji] = 0). Then the usual bosonization relation 

^ aj {x) = a -^ Fa . e -i(N aj -l/2)2nx/L ( 4 ) 

holds as operator identity, where the Klein factors |^] 
F aj = Vo^(0)e^C°J satisfy (see §) [F aj ,N a , r ] = 
S a a'Sjj'F a j, [F, (f>]=0 and {F aj ,Fl lj ,} = 26 aa >6 jj <. Thus 
F a j , i*^ ladder between the N a j , A^q.^ =F 1 Hilbert spaces 
without creating particle- hole excitations, and ensure 
proper i/i^tp' anticommutation relations. 

To exploit the conserved quantities in the Af y basis, we 
now use the transformation (^]) to define new bose fields 
(f y . Writing H in terms of these 



uqaj 
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[via (Q)], only ip x and ip s couple to the impurity Q: 
#o = A L £~A^ + ]T 9b\ y b qy , (5) 



y,n g >0 



H z = \ z A L S z Af s + X Z A L S Z ^ - b\ s ) , 



(6) 



n„>0 



H_l = ^ e^ >S_ V Pl-i^-e^ ) + H.c. (7) 
2o — ' 

To eliminate P z , make the EK || unitary transformation 
H' = UHU\ with U{\ z ) = e * A ^^(°). This yields 



H' h = H h , (H + H z )' = H + \ z A L Af s S z + const., S' ± = 
e ±i\zip s (o) an( j p s incurs a phase shift: 

Ucp a (x)U^ = cp g (x) - X z TtS z sga(x) ee ip s (x). (8) 

We henceforth focus on the EK line of fixed \ z = 1. 
Here tp s decouples from S± , and by (jj) and (||) the ip a j 's 
have phase shifts ±7r/4. Since this is just the value known 
for the NFL fixed point the A_L-induced crossover 

between the free and NFL fixed points can be studied on 
the EK line || by solving H' by refermionizing. 

Refermionization. — We first have to define Klein fac- 
tors for the Af y basis. Since an "off-diagonal" product 
F^jF a /ji acting on any state \N) just changes some of its 
N a j ( and hence Af y ) quantum numbers, we write 



(9) 



thereby defining new Klein factors T y , satisfying 

■t 



28,,, 



[T yi Ny'] = byy'Fy, [f , <f\ = and {T y: Ty,} 
Formally, these operators act on an extended Fock space 
of states with arbitrary Af E (Z + P/2) 4 . Its 
physical subspace contains only those states that obey 
(^), and by (0) it is closed under the pairwise action 
of TyS,. This simple construction for keeping track of J\f y 
quantum numbers is the main innovation of this Letter. 
Next we define a pseudofermion field ip x (x) § by 

i> x { X ) EE a -l/2^ e - J (A f .-l/2)2Wi e -^W j ( 10 ) 

and expand it as VAl J2 k e ~ %kXc kx> by analogy with (||) 
and (0), which imply {c hx ,c\, x } = 5 55 ,. In the c kaj ba- 
sis, the c^'s create highly non- linear combinations of 
electron-hole excitations, as is clear from their explicit 
definition, via ip x and T x , in terms of the Cfc a j's. Since 
A4 G Z + -j, we note that ip x has a P-dependent bound- 
ary condition, implying k = A^(nj. 

ther that A L (N 2 J2 + J2 q>0 n q b\ x b qx 



^-^), and fur- 

= H 0x + P/8, 

: means normal 
\ x c- kx : = J\f x - P/2. We 



and 



where H 0x = J2k^ 
ordering of c^'s, with J^k 
further define the "local pseudofermion" Cd = J-JS-, 
implying c^c^ = S z + ^. Eliminating J\f s in the sub- 
space S using A/" s = St + \ — c\cd we can rewrite H' 
as H csf (b c , bf,b s ,Af c ,Af f ) + H x + E G , where H csf has a 
trivial spectrum and H x is quadratic: 

H x - e d c\c d + H 0x + \fA L sY j {c\ x + c- kx ){c d - 4), 

k 

E G = A L [\(S 2 T -\) + P/8] - \K + h e (S T + |). 

Here T = A^/4a and e d = hi — h e is the spin flip en- 
ergy cost. As first noted by EK M, who derived H' 
for L — > oo, impurity properties show NFL behavior since 
"half the pseudofermion", (c d + c\), decouples. 
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Diagonalizing H x . — To study the NFL behavior of 
electron properties, caused by nonconservation of J\f x , 
we diagonalize H x . First define pseudofermions hav- 
ing all non-negative energies: aj, = -^=(0^ + c_ £fl; ) and 



1/2, which imply that Aj_ grows 



fa = -rm^kx ~ c_- h ) for k > 0; if P = 1 then a 



and ay = Q (or c\) for > (or < 0). Then the 
/%'s decouple in H x , and a Bogoliubov transformation 



E„=d,fc £r/=± BenMl + va n ) yields (1 



Ed 



I) + (/j!ft + 1), 



47rre/(e 2 -^) = - cot tt(s/A l - P/2) . (11) 



Equation for the pseudofermion eigenenergies £ im- 
plies that each k smoothly evolves into a corresponding 
e(k) as T is turned on. Since s(k) ~ fc + -g^ (or ~ fe) 
for fe <C (or ^>) r, we see very nicely that the spectrum's 
low- and high-energy parts are strongly and weakly per- 
turbed, respectively, with crossover scale T K ^ T §. 

As mentioned above, the pseudofermions act on an ex- 
tended Fock space. To identify which eigenstates \E) of 
H' are physical, note that each has to adiabatically de- 
velop, as r increases from 0, from some state obeying 
the free gluing condition (||). The latter can be shown 
pj to develop into the general gluing condition (GGC) 



mod2\E) must 



@ that (E\ [Ee>o4&e+Ek>0^fa\ 

be equal to [j\[ c + Nf - (St + 5 + f - Pd)]mod2, where 
P d = 0(1) for e d > (< 0). The GGC and Eqs. @ to- 
gether constitute an exact analytical solution of the 2CK 
model at the EK-line for arbitrary Aj_, hi and h e . 

Relation to RG methods. — Our exact solution allows 
us to implement Anderson's "poor man's scaling" and 
Wilson's NRG treatments of the Kondo problem analyti- 
cally, thus illustrating the main idea behind both, namely 
to try to uncover the low- energy physics via an RG trans- 
formation. In the first, the RG is generated by reducing 
(at fixed L, usually = 00) the bandwidth while adjusting 
the couplings to keep the dynamical properties invariant. 
Since the cut-off used when bosonizing is 1/a (~ pp) and 
a occurs in H' only through T, the scaling equations p| 
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FIG. 1. All eigenenergies AE = (E - E min )/A L < 1 
(degeneracies in brackets) of the full H as functions of (a) 
A z e [0,1] at T = e d = 0; (b) T / A L e [0, 00] at \ z = 1, 
£d = 0; (c) \e d \/T G [0, 3.5] at fixed V/Al > 1, A* = 1. 
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under rescaling |T|], are exact along the EK-line. Renor- 
malizing the spin-flip vertex, possible only approximately 
in the original Ckaj basis by summing selected diagrams, 
thus becomes trivial after bosonizing and refermionizing, 
which in effect resums all diagrams into a quadratic form. 

Wilson's NRG f§| is, in effect, a finite-size scal- 
ing method which increases (at fixed bandwidth and 
couplings) the system size, thus decreasing the mean 
level spacing and pushing ever more eigenenergies down 
into the spectrum's strongly-perturbed regime below Tk- 
Each RG step enlarges the system by order A > 1 by 
including an extra "onion-skin shell" of electrons, then 
rescales H — > AH to measure energy in units of the new 
reduced level spacing. We can mimick this by transform- 
ing L -f L' = AL (thus T/A L -> Ar/A L ) and plotting 
the spectrum in units of A' L = jr. 

Figure [lj displays (E — E rn i n )/ Al for the lowest few 
\E) that satisfy the GGC. Figure |](a) shows the evolution 
of the spectrum toward the EK-line for X z G [0, 1] at T = 
e d = [i.e. free fermions, phase-shifted by ±A 2 7r/2 in 
the spin sector, see (0)]. Figure [l](b) shows its fur- 
ther evolution on the EK-line for T/A^ G [0,oo] at 
X z = 1, e d = 0. Decreasing A^ at fixed T yields an 
NRG-like crossover spectrum that for A^ — > indeed re- 
produces the NRG's universal NFL fixed point spectrum 
(irrespective of the specific T value, illustrating the 
irrelevance of spin anisotropy |^]). This NFL spectrum 
also agrees with that found by AL using a so-called fusion 
hypothesis [0] , which our GGC thus proves simply and di- 
rectly (in contrast to the CFT proof of Ref. [14(b)]). Note 
that the ground state (with degeneracy 2) has entropy 
In 2, as it must for finite L jl5) (in contrast, the celebrated 
result \ In 2 requires taking L — ► 00 before T — * 0). 

Next we illustrate Wilson's program of extracting the 
most relevant operator's dimensions from the L depen- 
dence of the finite-size corrections, SE(L) = E(Y/ Al) — 
E(oo), to the universal NFL spectrum: For e d = 0, 

thus on the EK line the least 



Eq. (ITT]) gives SE 



j_ 

irrelevant operator has dimension 1; but perturbative 
corrections in A 2 1 ■ ,u - A 1 



1 yield 4£ 



thus the gen- 



eral leading irrelevant operators (absent on the EK line) 
have dimension ^ [p| JT^ , p^ . Next, turning on a local 
field e d — hi, we find from ( ]TT"| ) that for hi <C h c = 
yT/Z the NFL spectrum is only slightly affected, while 
for h c <C hi <§; r the spectrum has three distinct regions: 
It is Fermi-liquid-like 0] (with uniform level spacing) for 

e ^ hK = -fr and e ^ T, and NFL- like (nonuniform 
level spacings) for hx <C e <C T. Both the L depen- 
dence of h c and the hi dependence of the crossover scale 
hx show that the local magnetic field is relevant, with 
dimension — |; it causes a crossover, showsn in Fig. |l|(c), 
to a Fermi- liquid spectrum for all states with e <C hx- 

For r/Ai — > cxd, h — * 0, we find logarithmic diver- 
gences for the susceptibility, x ~ j^t^ u (^L): an d t ne 
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Af x fluctuations, (A/" 2 ) « 4* ln(rL) (with (Af x ) = 0). 
Both are clear signs of 2CK NFL physics: the first shows 
that no spin singlet is formed due to "overscreening" , the 
second how strongly this perturbs the electron sea. 

Relation to CFT. — Recent CFT § and scaling § ar- 
guments showed that the NFL regime can be described 
by free boson fields. This can be confirmed very easily 
by finding the scattering state operators cl [and field 

Tp x {x)] into which the free cj^'s [V4( x )] develop when 
r is turned on adiabatically as e vt T (at Ed — 0), and 
deducing from these the behavior of the (p y fields. In 
the continuum limit [L — > oo, then (A^ <c) r\ — > + ], 
the cl 's obey [ [l6| the Lippmann-Schwinger equation 



4 

k'l 



] = k 



st 



ir )^kx ~ C L)> wnicn g ives l^. 



dk'2Tk(ct+c_- klx ) 



[(k + in)(k + iAnT) - e 2 d ] (k - k' + if]) ' 



To find the asymptotic behavior (\x 
we may take k/T 




-> oo) of Wx{x) 
0: this gives 



dk'e 1 



[cl,6(x) - c_- k ,J(-x)] 



Adopting AL's notation of L- and R- movers, tp x (x) = 

0( x )$Il( x ) + (- x )i>l R ( x ), this S ives ~i>xR ~ #L ~ 
To translate this into "boundary conditions" on the 
<p y boson fields, we write 4> x l/r = F X L / 'roT 1 ! 2 '■e~ t v* 1 '/ 'R 
and note that tp c , iff decouple and (p s is phase-shifted 
as in (||). Thus the free and scattering boson fields are 
asymptotically related (with T) c ,T] s ,r}f = 1 = — rj x ) by 



(VytfyR - nS Z 5y S ) ~ (ifyL + irS z 5ys) 



9y 



(12) 



while r]y(J : yR ) r >v = T yL = T y for y = s,f,x. This cen- 
tral result, first found in Ref. (with different phases 
since Klein factors were neglected), shows that the NFL 
regime can be described by boson fields (p y L/R that are, 
asymptotically, free (with only a trivial S z dependence). 

Next we consider the 16 bilinear fermion currents 
j« A = :^ aj (T^ A ) aa > dj 4 a > f : (with T c 00 = ±55, Tf = 
\a a 5, T° A = ±5a A , T X A = \a a a A ), for which (|iT 



yields |Tj the "boundary conditions" ~ VyJyi- ^or 
y = c, s, /, these express the reemergence at the NFL 
fixed point of the full U(l) x SU(2) 2 x SU(2) 2 Kac- 
Moody symmetry assumed by AL; and for y = x they 
are just what AL derived using their fusion hypothe- 
sis. Since these boundary conditions fully determine all 
AL's CFT Green's functions (J], the boson approach will 
identically reproduce them too, if one proceeds as fol- 
lows: to evaluate (Vv,(l) . . . ^/^/(l'))) simply insert (||), 

rewrite the result in terms of <p v l/r and ^F v l/r, and 
combine ( |l2|) with standard free-boson results like 



/ e -i\<p yR (x) e i\' ip y , L (x')\ 



S ,L-!/2(»JyA-A') 2 

(13) 



All asymptotic NFL behavior of electron Green's func 



Hons arises from the fact that rj x = —1, combined with 
relations like (|l3|); it directly yields, e.g., the so-called 
"unitarity paradox" @ (^ajR{x)^ a , j, L {x')) ~ (for L -> 
oo, then \x' — x\ — * oo). Note, though, that probability is 
not lost during scattering: tpl (x) shows that each pseudo- 
particle c^_, incident from x > is "Andreev-scattered" , 
emerging at x < as pscudo/iofe c_y x , orthogonal to 
what was incident; this very NFL-like behavior dramati- 
cally illustrates the effects of N x nonconservation. 

To find AL's boundary operators in terms of the ifyS 
|^|,^lj, one calculates the operator product expansion of 
4>RajYLa'j'- Since rj x = —1, all terms contain a factor 
e ±vp y (y = S J or x) with dimension |; this ultimately 
causes the famous T 1 / 2 in the resistivity ^^,0. 

In conclusion, finite-size bosonization allows one (i) to 
mimick, in an exact way, the strategy of standard RG 
approaches, and (ii) to recover with remarkable ease all 
exact results known from CFT for the NFL fixed point. 
It thus constitutes a bridge between these theories. 
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